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Abstract 

In this paper, a new control scheme, called additive-state-decomposition-based tracking control, 
is proposed to solve the tracking (rejection) problem for rotational position of the TORA (a nonlinear 
nonminimum phase system). By the additive state decomposition, the tracking (rejection) task for 
the considered nonlinear system is decomposed into two independent subtasks: a tracking (rejection) 
subtask for a linear time invariant (LTI) system, leaving a stabilization subtask for a derived nonlinear 
system. By the decomposition, the proposed tracking control scheme avoids solving regulation 
equations and can tackle the tracking (rejection) problem in the presence of any external signal 
(except for the frequencies at ±1) generated by a marginally stable autonomous LTI system. To 
demonstrate the effectiveness, numerical simulation is given. 

Index Terms 

TORA, RTAC, Nonminimum phase, Additive state decomposition. 

I. Introduction 

The tracking (rejection) problem for a nonlinear benchmark system called translational 
oscillator with a rotational actuator (TORA) and also known as rotational-translational ac- 
tuator (RTAC) has received a considerable amount of attention these years Q]|-[[_]|. Some 
results were presented concerning the tracking (rejection) problem for general external sig- 
nals |ED,||3). However, the proposed control methods cannot achieve asymptotic disturbance 
rejection. Taking this into account, the nonlinear output regulation theory was applied to 
track (reject) external signals generated by an autonomous system. In this case, asymptotic 
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disturbance rejection can be achieved. By using different measurement, the tracking (rejection) 
problem for translational displacement of the TORA were investigated @]|-[[6]|. Readers can 
refer to [6] for details. Based on the same benchmark system, some other work was also 
presented concerning the tracking (rejection) problem for rotational position by nonlinear 
output regulation theory 0,|[S). For the two types of tracking (rejection) problems, regulator 
equations have to be solved and then the resulting solutions will be further used in the 
controller design. However, the difficulty of constructing and solving regulator equations 
will increase as the complexity of external signals increases. Moreover, it may fail to design 
a controller if regulator equations have no solutions. These are our major motivation. 

In this paper, the tracking (rejection) problem for rotational position of the TORA as [JVjl , [[HI 
is revisited by a new control scheme called additive-state-decomposition-based tracking 
control, which is based on the additive state decomposition^. The proposed additive state 
decomposition is a new decomposition manner different from the lower-order subsystem 
decomposition methods. Concretely, taking the system x (t) = f (t, i),i6M" for example, 
it is decomposed into two subsystems: x\ (t) = fi (t, Xi, x 2 ) and x 2 (t) = f 2 (t, x±, x 2 ), where 
x\ G M ni and x 2 G MJ 12 , respectively. The lower-order subsystem decomposition satisfies 

n = rii + n 2 and x = x\ © x 2 . 

By contrast, the proposed additive state decomposition satisfies 

n = rii = n 2 and x = x\ + x 2 . 

In our opinion, lower-order subsystem decomposition aims to reduce the complexity of 
the system itself, while the additive state decomposition emphasizes the reduction of the 
complexity of tasks for the system. 

By following the philosophy above, the original tracking (rejection) task is 'additively' 
decomposed into two independent subtasks, namely the tracking (rejection) subtask for a 
linear time invariant (LTI) system and the stabilization subtask for a derived nonlinear system. 
Since tracking (rejection) subtask only needs to be achieved on an LTI system, the complexity 
of external signals can be handled easier by the transfer function method. It is proved that 
the designed controller can tackle the tracking (rejection) problem for rotational position of 

'in this paper we have replaced the term "additive decomposition" in (9) with the more descriptive term "additive state 
decomposition". 
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the TORA in the presence of any external signal (except for the frequency at ±1) generated 
by a marginally stable autonomous LTI system. 

This paper is organized as follows. In Section 2, the problem is formulated and the 
additive state decomposition is recalled briefly first. In Section 3, an observer is proposed to 
compensate for nonlinearity; then the resulting system is 'additively' decomposed into two 
subsystems; sequently, controllers are designed for them. In Section 4, numerical simulation 
is given. Section 5 concludes this paper. 

II. Nonlinear Benchmark Problem and Additive State Decomposition 

A. Nonlinear Benchmark Problem 

As shown in Fig.l, the TORA system consists of a cart attached to a wall with a spring. 
The cart is affected by a disturbance force F. An unbalanced point mass rotates around the 
axis in the center of the cart, which is actuated by a control torque N. The translational 
displacement of the cart is denoted by x c and the rotational position of the unbalanced point 
mass is denoted by 9. 




Fig. 1. TORA system configuration 1 1 1 



For simplicity, after normalization and transformation, the TORA system is described by 
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the following state-space representation |0Q: 

x\ = X2 (la) 
x 2 = —X\ + e sin x 3 + F d (lb) 



±3 = £4 (lc) 

£ COS X 3 

1 — e 2 cos 2 £3 



c wjo j/3 . . 

x 4 = m-- 5 5 — F d ,x(0) = xo (Id) 



where < £ < 1, 2 = [ xi x 2 x 3 x 4 ] T G M 4 , x 3 = 6>, F d G R is the unknown 
dimensionless disturbance, u G R is the dimensionless control torque. In this paper, the 
tracking (rejection) problem for rotational position of the TORA as ||7),|[E) is revisited. 
Concretely, for system (OQ), it is to design a controller u such that the output y (t) = x 3 (t) — > r 
as t — > 00, meanwhile keeping the other states bounded, where r G (— vr/2 , 7r/2) is a 
known constant. Obviously, this is a nonlinear nonminimum phase tracking problem, or 
say a nonlinear weakly minimum phase tracking problem. For system ©, the following 
assumptions are imposed. 
Assumption 1. The state x can be obtained. 

Assumption 2. The disturbance F d G R is generated by an autonomous LTI system 

«; = ^,F d = C> (2) 

where S = —S T G ]R mxm , G M m are constant matrix, wel™, and the pair (Cj, S) is 
observable. 

Remark 1. If all eigenvalues of S have zero real part, then, in suitable coordinates, the 

matrix S can always be written to be a skew-symmetric matrix. The matrix S in previous 

u 

literature on the output regulation problem is often chosen in a simple form S = 

-u 

where u is a positive real H-flH. In such a case, F d is in the form as sin(±ait) and the 
solution to the regulator equation is easier to obtain. However, this is a difficulty when S is 
complicated. 

B. Additive State Decomposition 

In order to make the paper self-contained, the additive state decomposition is recalled 
here briefly. Consider the following 'original' system: 

f(t,x,x) = 0,x(0) = x (3) 
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where x E W 1 . We first bring in a 'primary' system having the same dimension as ©, 
according to: 

fp(t,Xp,x p ) = 0,x p (0) = x pfi (4) 

where x p E M n . From the original system © and the primary system © we derive the 
following 'secondary' system: 

/ (t, x, x) - f p (t, x p , x p ) =0,x (0) = x (5) 

where x p E IR n is given by the primary system ©. Define a new variable x s E W 1 as follows: 

x g — x x>p> 

Then the secondary system © can be further written as follows: 

/ (tj is &pi *£s %p) fp \b-> %pi 0, x s (0) Xq 3-p,o- (7) 

From the definition ©, we have 

x (i) = xp (t) +x,(t),t> 0. (8) 

Remark 2. By the additive state decomposition, the system © is decomposed into two 
subsystems with the same dimension as the original system. In this sense our decomposition 
is "additive". In addition, this decomposition is with respect to state. So, we call it "additive 
state decomposition". 

As a special case of ©, a class of differential dynamic systems is considered as follows: 

X = f (t,x) ,X (0) = X , 

y = h(t,x) (9) 

where x E W 1 and y E R m . Two systems, denoted by the primary system and (derived) 
secondary system respectively, are defined as follows: 

Xp fp (t, Xp) , Xp (0) Xp Q 

y p = h p (t,x p ) (10) 

and 

X s f (t, Xp -\- X s ) fp (t : Xp) , X s (0) Xq ^p,0) 

y s = h(t,x p + x s ) - hp(t,x p ) (11) 
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where x s = x — x p and y s = y — y p . The secondary system (fTTT) is determined by the original 
system © and the primary system (flOl) . From the definition, we have 

x (t) = x p (t) + x s (t) ,y(t)= y p (t) +y s (t),t> 0. (12) 

III. Additive-State-Decomposition-Based Tracking Control 

In this section, in order to decrease nonlinearity, an observer is proposed to compensate for 
the nonlinear term ^2^2 x Fa,- After the compensation, the resulting nonlinear nonminimum 
phase tracking system is decomposed into two systems by the additive state decomposition: 
an LTI system including all external signals as the primary system, leaving the secondary 
system with a zero equilibrium point. Therefore, the tracking problem for the original system 
is correspondingly decomposed into two subproblems by the additive state decomposition: a 
tracking problem for the LTI 'primary' system and a stabilization problem for the secondary 
system. Obviously, the two subproblems are easier than the original one. Therefore, the 
original tracking problem is simplified. The structure of the closed-loop system is shown in 
Fig.2. 



Tracking Controller for 
Primary system 



Stabilized Controller 
for Secondary System 



+ 

v. 



x p ,x s 



Decomposed System 
Observer 



W- 



TORA 



Nonlinearity 
Compensation 



Fig. 2. Structure of the closed-loop system 



A. Nonlinearity Compensation 

First, in order to estimate the term ,, £ o OSX i — F rf , an observer is designed, which is stated 

' 1— £^ COS^ £3 O ' 

in Theorem 1. 
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Theorem 1. Under Assumptions 1-2, for system (OQ), let the observer be designed as follows 

e cos x 3 



1 — e 2 cos 2 x 3 



w = Sw + l\ 
x 4 - 

F d = kCjw, w (0) = 0, x 4 (0) = 



/ 2 (X4 - X 4 ) - <!- C rf 10 + M 



1 — e 2 cos 2 x 3 



(13) 



where l\,l 2 > 0. Then lim T 
Proof. See Appendix A. □ 



£ COS 3:3 



: F d (t) = 0, where F d ^F d - F d . 



By using the observer (PT3T) . the controller it in (OQ) is designed as follows 



u = K x + v + 



e cos £3 



1 — e 2 cos 2 x 3 



where K G IR 4 and »Gl will be specified later. Then the system CD becomes 

Xi = x 2 

x 2 = —xi + e sin x 3 + F d 
x 3 = £4 



£4 = K x + v + 



e cos x 3 



1 — e 2 cos 2 £3 



F d) x(0) = £ - 



(14) 



5. Additive State Decomposition of Original System 

Introduce a zero term eD (C + a5) x — e (y + ay) = into the system (|14l) . where a > 0, 
-B = [ 1 ] T , C = [ 1 ] T and (C + glB) t x = y + ay. Then the system 
(fl"4l) becomes 

i = Ax + Bv + <j> (y, y) + DF d + ip 
y = C T x, x (0) = xq 



(15) 



where 



D 






1 





-1 














1 











" 






1 








A 


(y,y) 

















A = A + BK t + eD(C + aB) J 





e sin y — e (y + ay ) 











£ COS X3 jp 

1 — e 2 cos 2 Z3 <* 



(16) 
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The additive state decomposition is ready to apply to the system (fl"5T) . for which the primary 
system is chosen to be an LTI system including all external signals as follows 

x p = Ax p + Bv p + d + cp 

y p = C T x p , x p (0) = x Q (17) 

where d = (j) (r, 0) + DF d . Then, according to the rule (PTTI) . the secondary system is derived 
from the original system (fl"5l) and the primary system (PT71) as follows 

x s = Ax s + Bv s + 4> (y p + y s , y p + y s ) - (r, 0) 
y s = C T x s ,x s (0) = (18) 
where v s = v — v p . According to (PT2l) . we have 

x = x p + x s and y = y p + y s . (19) 

Remark 3. The pair (A Q ,B) is uncontrollable, while the pair (A + eD{C + aBy ,B) 
is controllable. Therefore, there always exists a vector K such that A = A + BK T + 
eD (C + aB) T is a stable matrix. 

Remark 4. If y p = r and y p = 0, then (x s ,w s ) = is a zero equilibrium point of the 
secondary system (PT8l . 

So far, the nonlinear nonminimum phase tracking system (PT5l) is decomposed into two 
systems by the additive state decomposition, where the external signal d + ip is shifted to 
(fTTT) and the nonlinear term (■) is shifted to (PT81) . The strategy here is to assign the tracking 
(rejection) task to the primary system (fT71) and stabilization task to the secondary system 
(fT8l . More concretely, in (fT71) design v p to track r, and design w s to stabilize (fT8l . If so, 
by the relationship (TT9l , y can track r. In the following, controllers v p and w s are designed 
separately. 

C. Tracking Controller Design for Primary System 

Before proceeding further, we have the following preliminary result. 
Consider the following linear system 

zi = S z zi + A 12 e z 

z 2 = A%\Z\ + A 22 z 2 + di + ipi 

e z = C^z 2 + d 2 + <p 2 ,z{0) = z (20) 
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where S* G 



ami xmi 



is a marginally stable matrix, A\ 2 G R mi , C e G 



P'"2 



A 2 i g 



pm2 xrti2 



^22 g 



Dm2 xmi 



, «i G 



omi 



2^,di,^i G 



P' m 2 



Dm 1+ m 2 and e ^d 2 ,(f 2 £ 



Lemma 1. Suppose i) (pi(t) is bounded on [0, 00) and lim = 0, i = 1,2, ii) 

every element of di (t) , d 2 (t) are bounded on [0, 00) and can be generated by w z = S z w z , 
d z = Cjw z with appropriate initial values, where C z G M mi , iii) the parameters in (1201) 
satisfy 



max Re A (A z ) < 0, A z 



A 



21 



.4 



22 



(21) 



Then in (1201) lime z (£) = 0, meanwhile keeping z x (£) and z 2 (£) bounded. 

t— ¥00 

Proof. See Appendix B. □ 

Define a filtered tracking error to be 

e P = y P + ay p = (C + oB) x p - r 



(22) 



where y p = y p — r, r = and a > 0. Let us consider the tracking problem for the primary 
system (fT7l) . With Lemma 1 in hand, the design of v p is stated in Theorem 2. 

Theorem 2. For the primary system (PT71) . let the controller v p be designed as follows 



(23) 



max Re A (A a ) < 0, A a 



(24) 



£ = S a £ + L x e p 
v v (?) x pi r ) = L^Xp + L3 £ 
where S a =diag(0, S) , L x e R m+1 , L 2 G M 4 and L 3 G M m+1 satisfy 

5 a MC + aBf 

Then lim y p (£) = r and lim y p (£) = meanwhile keeping x p (£) and £ (£) bounded. 

t— ¥00 t— >oo 

Proof. Incorporating the controller (1231) into the primary system (fTTT) results in 

I = ^ + L x e p 

x p = (A + BL\) x p + BLl^ + d + cp 
e p = (C + aB) T x p — r 

where the definition (l22l) is utilized. Moreover, every element of d and r can be generated by 
an autonomous system in the form w a = S a w a , d a = C^w a with appropriate initial values, 
where C a — [ 1 Cj ] T . By Lemma 1, if (|24l) holds, then lime„ (£) = meanwhile keeping 
Xp (£) and £ (£) bounded. It is easy to see from (l22l) that both y p and y p can be viewed as 
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outputs of a stable system with e p as input. This means that y p and y p are bounded if e p is 
bounded. In addition, limiL (t) = and limy„ (£) = 0. □ 

In most of cases, the controller parameters Li,L 2 and L 3 in (1231) can be always found. 
This is shown in the following proposition. 

Proposition 1. For any S = —S T without eigenvalues ±j, the parameters 

L x = [\ cJ] T ,L 2 = --C-B--H--K,L 3 = --L 1 (25) 

a a a a 

can always make max Re A (A a ) < 0, where H = [ e 1 ] T - 
Proof. See Appendix C. □ 

Remark 5. Proposition 1 in fact implies that, in the presence of any external signal 
(except for the frequencies at ±1), the controller (1231) with parameters (1251) can always make 
lim y p [t) = r and lim y p (t) = meanwhile keeping x p (t) and £ (£) bounded. In other words, 
the disturbance like sint cannot be dealt with, which is consistent with Q. If the external 
signal contains the component with frequencies at ±1, then such a frequency component can 
be chosen not to compensate for, i.e., S a in (1231) will not contain eigenvalues ±j. 

D. Stabilized Controller Design for Secondary System 

So far, we have designed the tracking controller for the primary system (fTTT) . In this 
section, we are going to design the stabilized controller for the secondary system (PT8T) . It can 
be rewritten as 

&l,s — x 2,s 

£2,8 = -Xi,s + £ sin (x 3)S + r) - e sin r + g 

2-3, s ^4,s 

X4,s = K T x s + v s ,x s (0) = (26) 

where g = e sin (y p + x 3:S ) — e sin (r + x 3 s ) — e (y p + ay p — r) . Our constructive procedure 
has been inspired by the design in [3J. We will start the controller design procedure from 
the marginally stable (xi >s , X2, s )-subsystem. 

Step 1. Consider the (xi )S , x 2iS )-subsystem of (1261) with x 3jS viewed as the virtual control 
input. Differentiating the quadratic function V\ = \ (x\ s + x\ s ) results in 

Vi = £X2, S [sin (x 3iS + r) - sin r] + ex 2 , s g- (27) 
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Guided by the state-feedback design ifTOl . we introduce the following "Certainty Equivalence" 
(CE) based virtual controller 

x 3 , s = -ba.ta.nx 2 , s + x 3 >a . (28) 

Then 



where 



&2,s 



X 2r 



^ . . -6atanx 2 «\ /-6atanx 2 . + 2r\ , 
-xi iS + 2e sin ( ) cos ( ) + g 



g' = e sin (r — fcatanx 2s + x' 3 J — e sin (r — f>atanx 2s ) + g. 



(29) 



(30) 



In order to ensure cos ^ batanx ^ a+2r j > q ; (j-^ parameter 6 is chosen to satisfy < b < 
2 (1 — 2 |r| /it) . Since r G (— it /2 , 7r /2) is a constant, 6 always exists. The term CE is used 
here because x' 3s = in (1281) makes V\ in (1271) negative semidefinite as g = 0. 

Step 2. We will apply backstepping to the (x' 3 s , x^ s ) -subsystem and design a nonlinear 
controller v s to drive x' 3 s to the origin. By the definition (1281) . Xg = x 3jS + 6atanx 2jS . Then 
the time derivative of the new variable x' 3 s is 



x' 



3..s 



x 4 s + ^ + b 



1 + 4/ 



(31) 



where ^ = [— xi, s + £ sin (x 3)S + r) — e sin r] . Define a new variable x 4s as follows 



x' 



4,s 



4,8 + X i,s + i>- 



(32) 



Then (1311) becomes 



ar 



3,8 



4,s 



1 + ^9 a 



By the definition (1321) . the time derivative of the new variable x 4 s is 



where 



1 



-x(, , +x' As + b x g + fT T x s + v s + ijj, 
1 + 



-26x 2 



[— ^i,s + £ sin (x 3jS + r) — e sin r] 



+ 6 



1 + x 2 



[—^2,s + £ cos (x 3jS + r) x 4 ,. 



2,s 
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Design v s for the secondary system (|26l) as follows 

v s (x p , x s , r) = x' 3s - 2x' 4yS - K T x s - tp. (33) 
Then the (x' 3 s , x' 4 s ) -subsystem becomes 

X 3,s — ~ X 3,s + X 4,s + ^^_j_ x 2 9 



1 

''V 



#m = -< )S + h YTx r9 - (34) 



It is easy to see that lim £3 (t) = and lim £4 (t) = as lim (7 (t) = 0. 

t— >oo ' t— s-oo ' t— >oo 

We are now ready to state the theorem for the secondary system. 

Theorem 3. Suppose lim y p (t) = r and lim ?/„(£) = 0. Let the controller v s for the sec- 

t— s-oo t— s-oo 

ondary system (1261) be designed as (1331) . where < 6 < 2 (1 — 2 |r| /7r). Then lim ||x s (t)\\ = 

t— >OD 

meanwhile keeping x s (t) bounded. 
Proof. See Appendix D. □ 

E. Controller Synthesis for Original System 

It should be noticed that the controller design above is based on the condition that x p and 
x s are known as priori. A problem arises that the states x p and x s cannot be measured directly 
except for x = x p + x s . By taking this into account, the following observer is proposed to 
estimate the states x p and x s , which is stated in Theorem 4. 

Theorem 4. Let the observer be designed as follows 

x s = Ax s + Bv s + <p (y, y) - <p (r, 0) 

x p = x — x s , x s (0) = (35) 

where A is stable. Then x p = x p and x s = x s . 

Proof. Since x = x p + x s , we have y — y p + y s . Consequently, (1331) can be rewritten as 

x s = Ax s + Bv s + 4>(y p + y s , y p + y s ) - (f) (r, 0) 

x p = x — x s , x s (0) = 0. (36) 

Subtracting (TT8T) from (l36l) results in 

i s = Ax s ,x s {0) = (37) 

where Then x s = x s . Furthermore, with the aid of the relationship 

we have x p = x p . □ 
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Remark 6. Unlike traditional observers, the proposed observer can estimate the states of the 
primary system and the secondary system directly rather than asymptotically or exponentially. 
This can be explained that, although the initial value x is unknown, the initial value of either 
the primary system or the secondary system can be specified exactly, leaving an unknown 
initial value for the other system. The measurement x and parameters may be inaccurate. In 
this case, it is expected that small uncertainties lead to x p close to x p (or x s close to x s ). 
From (1371) . a stable matrix A can ensure a small x s in the presence of small uncertainties. 

Theorem 5. Suppose that the conditions of Theorems 1-4 hold. Let the controller u in the 
system (OQ) be designed as follows 

i = S a i + L l [{C + aB) T x p -r] 

u = K x + v p (£, x p , r) + v s (x p , x s , r) + - - — F d (38) 

1 — e 2 cos 2 x 3 

where F d is given by (fl"3l) . x p and x s are given by (1331) . v p (•) is defined in (1231) . and v s (•) 
is defined in (1331) . Then limy (£) = r meanwhile keeping x and £ bounded. 

t— >oo 

Proof. Note that the original system (OQ), the primary system ([T71) and the secondary system 
(TTSl have the relationship: x = x p + x s and y = y p + y s . With the controller (1381) . for the 
primary system (TT7T) . lim?/p (t) = r meanwhile keeping x p and £ bounded by Theorem 2. On 

t— >oo 

the other hand, for the secondary system (TTSl . we have \imx s (t) = meanwhile keeping x s 

t— >oo 

bounded on [0, oo) by Theorem 3. In addition, Theorem 4 ensures that x p = x p and x s = x s . 
Therefore, lim y{t)=r meanwhile keeping x and £ bounded. □ 

t— >oo 

IV. Numerical Simulation 

In the simulation, set e = 0.2 and the initial value x = [0 0] T hi(OQ). The 
unknown dimensionless disturbance Fd is generated by an autonomous LTI system © with 
the parameters as follows 

,C d =[l f,w(0) = [ 0.02 f. 

The objective here is to design a controller u such that the output y (t) = x 3 (t) — > r = 0.5 
as t — > oo meanwhile keeping the other states bounded. 

The parameters of the observer (fT3l) are chosen as h = h = 10. In (fT6l) . the parameters of A 
are chosen as a = 1 and K = [ — £ —1 — 2 ] T - Then maxRe(A (A)) = —0.01 < 0. Since 
matrix S 1 does not possess the eigenvalues ±j, the parameters of the tracking controller (|23|) 
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of the primary system can be chosen according to Proposition 1 that L\ — [\ Cj } T , L 2 = 
4x i and L 3 = —L 1 . These make A a in (1241) satisfies maxReA(y4 a ) = —0.0084 < 0. The 
parameter b of the stabilized controller (1331) is chosen as 6 = 1.5(1 — 1 /ir) < 2 (1 — 2\r \ /it) . 

The TORA system (OQ) is driven by the controller (1381) with the parameters above. The 
evolutions of all states of (OQ) are shown in Fig. 3. As shown, the proposed controller u drives 
the output y (t) = rr 3 (t) — > 0.5 as t — > oo, meanwhile keeping the other states bounded. 



0.2 




20 40 60 80 100 120 140 160 180 




100 120 140 160 180 



Fig. 3. Evolutions of all states 

Unlike the output regulation theory, the proposed method does not require the regulator 
equations. If the disturbance Fd consists of more frequency components, i.e., S is more 
complicated, the designed controller above does not need to be changed except for the 
corresponding S and Cd- This demonstrates the effectiveness of the proposed control method. 
For example, we consider that the unknown dimensionless disturbance Fd is generated by an 
autonomous LTI system © with the parameters as follows 

1.5 
-1.5 



S = dmg (S 1 ,S 2 ),S 1 



2 






,s 2 = 


-2 





C d =[l 1 0] ,w(0) = [0 0.02 0.02]'. 

The controller in the first simulation is still applied to this case except for replacing S and 
C d (the dimension is changed correspondingly). Driven by the new controller, the evolutions 
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of all states of © are shown in Fig.4. As shown, the proposed controller u drives the output 
y (t) = x 3 (t) — > 0.5 as t — >• 00, meanwhile keeping the other states bounded. 
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Fig. 4. Evolutions of all states when disturbance is complicated 

V. Conclusions 

In this paper, the tracking (rejection) problem for rotational position of the TORA was 
discussed. Our main contribution lies in the presentation of a new decomposition scheme, 
named additive state decomposition, which not only simplifies the controller design but 
also increases flexibility of the controller design. By the additive state decomposition, the 
considered system was decomposed into two subsystems in charge of two independent 
subtasks respectively: an LTI system in charge of a tracking (rejection) subtask, leaving 
a nonlinear system in charge of a stabilization subtask. Based on the decomposition, the 
subcontrollers corresponding to two subsystems were designed separately, which increased 
the flexibility of design. The tracking (rejection) controller was designed by the transfer 
function method, while the stabilized controller was designed by the backstepping method. 
This numerical simulation has shown that the designed controller can achieve the objective, 
moreover, can be changed flexibly according to the model of external signals. 
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VI. Appendix 

A. Proof of Theorem 1 

The disturbance F d e R is generated by an autonomous LTI system © with an initial 
value w (0) . It can also be generated by the following system 

w = Sw,F d = hCjw (39) 

with the initial value j^w (0) . Subtracting (fldl) and (l39l) from (fT3l) results in 

w = Sw + k- - — C d x 4 

1 — e z cos z X3 

x 4 = -l 2 x 4 - h- - — C d w (40) 

1 — e z cos z xs 

where l\,l 2 > 0, and w = w — w. Design a Lyapunov function as follows 

Vi = -w w + -x 4 . 



Taking the derivative of V\ along (1401) results in 



V X = \w T {S + S T )w + hw T 1 *T X \ ( '/•' ^ 



e 2 cos 2 X3 



2 ecosx 3 T _ 

- - hx 4 - = — C d w. 

1 — e z cos z X3 

By Assumption 2, S + S T = 0. Then the derivative of V\ becomes 



Vi < -hx\ < 0. 

Since l 2 > 0, from the inequality above, it can be concluded by LaSalle's invariance principle 



021 that limx 4 (t) = and lim . T s ff Cjw (t) = 0. □ 

B. Proof of Lemma 1 

Before proving Lemma 1, we need the following preliminary result. 

Lemma 2. If the pair (A z , B z ) is controllable, then there exists a Co G IR m such that 

C T {sI m - A Z )^B Z = det { J m _ Az) 

where A z eR mxm and B z G W 11 . 
Proof. First, we have 

(sI m -A z )- 1 B z = N[s n ~ 1 ••■ 1 f /det (sl m -A z ) 

where N e R mxm . If the pair (A Z ,B X ) is controllable, the matrix N is full rank HB- We 
can complete this proof by choosing Cq = (iY _1 ) T [ • ■ • 1 ] T . □ 
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With Lemma 2 in hand, we are ready to prove Lemma 1. 

i) For the system (1201) . we have 

z (t) = e-^zo + [ e~ A ^- r) (d a + <p a ) (r) dr, t > 
Jo 

where d a = [ d%A7 2 dj } T and y? a = [ (p£ A\ 2 (pi ] T ■ Based on the equation above, since 
A (A z ) < and \\d a (t)\\, \[(p a are bounded on [0, oo), it is easy to see that \\zi (t)\\ and 
\\z 2 (t)\\ are bounded on [0, oo). 

ii) For the system (l20l) . the Laplace transformation of z (s) is 

Z (s) = (si - A z y X [d a (s) + (fa (s) + Zq] . 

Then z\ (s) = C a (si - A z )^ [d a (s) + ip a (s) + z ) , where C a = [ I mi miXm2 ] T - The 
condition A (A z ) < implies that the pair (S z ,Ai 2 ) is controllable. Otherwise, for the 
autonomous system z = A z z, the variable Z\ cannot converge to zero as S z is a marginally 
stable matrix. This contradicts with the condition A (A z ) < 0. Then by Lemma 1, there exists 
a C G W mi such that 

°« Zl (S) = det(s 1 I-S z ) 6z (S) • 

Then e z (s) can be written as 

e z (s) = det (si - S z ) C^z 1 (s) 

= Q (s) det (si - S z ) [d a (s) + (p a (s) + z ] ■ 

where Q (s) = C^C a (si - A z )~ l D a . Since every element of d a can be generated by w z = 

S z w z , d z = C T z w z , we have d a (s) = [Cj (si - 67 z ) _1 w Zfi (0)]( mi+TO2 ) x i, where w Zji (0) G R. 
Since (si — S z )~ x = det( -J_ g ) adj(g/ — S z ), e z (s) is further represented as 

e z (s) = Q (s) Cjw z (0) [Cjadj (si - S.) w z>i (0)] (tBl+ma)xl 

+ (si - A z y l [<p a (s) + z ] . (41) 



Since A (A z ) < and the order of A z is higher than that of S Z1 moreover \\tp a (i)\\ is bounded 
on [0, oc 
6D- □ 



on [0, oo) and lim \\(p a (t)\\ = 0, for any initial value w Z) i (0) , we have lime z (t) = from 

t— >oo " " t— s-oo 
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C. Proof of Proposition 1 

If we can prove that the following system 



i = Sat + Li (C + aB) T x P 



(A + BLl) x p + BLli 



(42) 



is asymptotic stable, then Re A (A a ) < holds. Choose a Lyapunov function as follows 



V 



r 2 



Xo „ + -p 



where p = x 3iP + ax^ p = (C + aB) x p . With the parameters L\ — [ 1 Cj } T , L 2 = 
-\C - B - \H - K and L 3 = -\L U the derivative of V along 02) is 

V = -p 2 . 

Define S = |x| V (x) = oj ; where x = [ £ T ] T - The remaindering work is to prove 
S = {x\x = 0}. If so, by LaSalle's invariance principle lfT2ll . we have lim = 

t— >oo 

0. Therefore, the system (l42l) with the parameters is globally asymptotically stable. Then 

ReA(y4 a ) < 0. 

Since V = =>• p = and a > 0, we have 5 = { = 0} . Let x be a solution 

belonging to 5 identically. Then, from (1421) . we have 







(43) 
(44) 
(45) 



From (1441. it holds that 



1]? 



1 






•) 


-1 





On the other hand, from (g3) and (1431) . it holds that 

1 



%2,p G ^2 



where matrix 5^ does not possess eigenvalues ±j as matrix S does not. Therefore X2, p G 
«Si Pi 5 2 = {0} and then ,t 1jP = 0, namely S = {x\ x p = 0} . 

Let x be a solution that belongs identically to S. Then Lj£ = 0. Since the pair (Cj, S 1 ) 
is observable, by the definition L\ = — i[ 1 ] T , the pair (L^, 5 a ) is observable as well. 
Consequently, we can conclude that £ = 0, namely 5=|x|x=[^ T ] T = j . □ 
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D. Proof of Theorem 3 

This proof is composed of three parts. 

Part 1. limx 3 (t) = 0, limx 4 (t) — and limg' (£) = as lim (7 (£) = 0. If lim y p (i) = 

t— >00 ' t— >O0 ' t—^OO t^rOO t— >00 

r and \imy„ (t) =0, then from the definition of 4>(y,y) , we have lim g [t] =0 no matter 

t— >oo t— >oo 

what y s is. According to this, it is easy from (l34l) to see that lim £3 a (i) = and lim £4 s (£) = 
when the controller v s for the secondary system d26l) is designed as (l33l . Then, in (T29l) . 
lim#' (t) = 0. 

t— >oo 

Par? 2. lima;i iS (t) = and limx2, s (t) = 0. Since < b < 2 (1 — 2 |r| /7r) , the derivative 
V\ in (1271) negative semidefinite when o' (t) = 0, namely, 

-foatana:2 S \ / — 6atanx 2 s + 2r\ 
V\ = 2ex 2 s sin cos — J < 0, 



2 / V 2 

where the equality holds at some time instant t > if and only if X2, s (t) = 0. By LaSalle's 
invariance principle [12J that limxi s (t) = and limx 2 . s (t) = when g' (t) = 0. Because 
of the particular structure of (xi )S , X2,<))-subsystem d29l) . by using |[T3l Lemma 3.6], one can 
show that any globally asymptotically stabilizing feedback for (xi )S , X2 lS ) -subsystem d29l) 
when g' (t) = achieves global asymptotic stability of (xi jS , a; 2 ,s) -subsystem (1291 ) when 
limg' (t) = 0. Therefore, based on Part 1, limxi )S (t) = and limx2, s (t) = 0. 
Par? i. Combining the two parts above, we have lim \\x s (t)\\ = 0. For the (xi 8 ,X2, 8 )- 

t—>oo 

subsystem and (x' 3s , x' 4 s )-subsystem, ||a; 8 (t) \\ is bounded in any finite time. With the obtained 
result lim ||x a (t)\\ = 0, we have \\x s (t)\\ is bounded on [0, 00). □ 

t— >oo 
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